SYMMETRIZED CHEBYSHEV POLYNOMIALS 
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Abstract. We define a class of multivariate Laurent polynomials 
closely related to Chebyshev polynomials, and prove the simple 
but somewhat surprising (in view of the fact that the signs of the 
coefficients of the Chebyshev polynomials themselves alternate) 
result that their coefficients are non-negative. We further show 
that a Central Limit Theorem holds for our polynomials. 



Introduction 

Let T„(x) and Um{x) be the Chebyshev polynomials of the first and 
second kinds, respectively. We define the following Laurent polynomi- 
als (which, for lack of a better name, we call the symmetrized Chebyshev 
polynomials: 
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Sni^Cj X\ , . . . , X^, 
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The function i?„ arises in the enumeration of conjugacy classes in the 
free group on k generators, more specifically, there is the following 
result: 



Theorem 0.1 ( |R,ivin99] ) . The number of cyclically reduced words of 
length k in homologous to ei[ai]-|-- ■ ■-|-er[ar] is equal to the coefficient 
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These results first appeared in the author's 1998 preprint "Growth in free groups 
(and other stories)," but seem to be of independent interest. The positivity result 
was preprint math.CA/0301210, but there appears to be no reason to separate it 
from the limiting distribution result, and many reasons to keep them together. 
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of xl^ ■ ■ ■ x^'' in 



(3) 2{V2^)'Rk{ 



r 



;xi,...,x,) + (r-l)[l + (-l)'=] 



Remark. The rescaled Chebyshev polynomial Tk{ax)/a'' is called the 
k-th Dickson polynomial Tk{x,a) (see |Schur73j ). 

The coefficients of generating functions of combinatorial objects are 
non-negative, so one is led to wonder for which values of c are the 
coefficients of Rn{c; xi, . . . , Xk) and 5'„(c; xi, . . . , Xk) nonnegative. In 
this note we give an essentially complete answer (see Theorems 12.11 
and 123) • We also write down an explicit formula (Eq. (fTH|) ) for the 
coefficients of i?„ and Sn- Furthermore, we use the positivity and tools 
of probability theory to analyze the distribution of the coefficients of 
Rn to prove a central limit theorem - Theorem 13.11 

1. Some facts about Chebyshev polynomials 

The literature on Chebyshev polynomials is enormous; IRivlinQO] is 
a good to start. Here, we shall supply the barest essentials in an effort 
to keep this paper self-contained. 

There are a number of ways to define Chebyshev polynomials (almost 
as many as there are of spelling their inventor's name). A standard 
definition of the Chebyshev polynomial of the first kind Tn{x) is: 

(4) Tn{x) = cos 72 arccos X. 

In particular, Tq{x) = 1, Ti(x) = x. Using the identity 

(5) cos(x + y) + cos(x — y) = 2 cos x cos y 

we immediately find the three-term recurrence for Chebyshev polyno- 
mials: 



The definition of Eq. (jH) can be used to give a "closed form" used in 
Section 121 



(6) 



T„+i(x) = 2xT„(x) - r„_i(x). 




Indeed, let x 



cos 6. Then 




while 
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SO that 

\ / \ n / \ n 

- — Vx^ — Ij + yx + Vx"^ — Ij = 3? exp(m^) = cos n9. 

We also define Chebyshev polynomials of the second kind Un{x), 
which can again be defined in a number of ways, one of which is: 

(8) Unix) = -1^t;+i(x). 

A simple manipulation shows that if we set x = cos 6, as before, then 

(9) c/„(.) = 5i!i(!i±i)^. 

smt^ 

In some ways, Schur's notation W„ = f/„-i is preferable. In any case, 
we have Uo{x) = 1, Ui{x) = 2x, and otherwise the f/„ satisfy the same 
recurrence as the T„, to wit, 

(10) Un+l{x) = 2xUn{x) - Un-l{x). 

From the recurrences, it is clear that for f = T,U, fn{~x) = (— 
or, in other words, every second coefficient of Tn{x) and Un{x) vanishes. 
The remaining coefficients alternate in sign; here is the explicit formula 
for the coefficient c^n-2m of a;"~^™ of T„(x) : 



i"i^ = (-ir^r ")2"— , m = 0,l,.... 



n 

n — m \ m / 12 



This can be proved easily using Eq. 

2. Analysis of the functions Rn and Sn- 

In view of the alternation of the coefficients, the appearance of the 
Chebyshev polynomials as generating functions in Section ?? seems 
a bit surprising, since combinatorial generating functions have non- 
negative coefficients. Below we state and prove a generalization. Re- 
markably, Theorems 12.11 and 12.21 do not seem to have been previously 
noted. 

Theorem 2.1. Let c > 1. Then all the coefficients of Rn{c] x) are non- 
negative. Indeed the coefficients ofx'^, x"""^, . . . , x"""*"^, x~"' are positive, 
while the other coefficients are zero. The same is true of Sn in place of 

Rn- 

Proof. Let be the coefficient of x^ in f/„((c/2)(x+ 1/x)). The recur- 
rence gives the following recurrence for the : 



(12) a:;+, = c(a^i + a'=+i^ 

Now we shall show that the following always holds: 
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(a) : > (inequality being strict if and only if n — /c is even). 

(b) : > max(a^I^, a^l^), the inequality strict, again, if and 
only if n — is even. 

(c) : > a^_2 (strictness as above). 

The proof proceeds routinely by induction; first the induction step 
(we assume throughout that n — A; is even; all the quantities involved 
are obviously otherwise): 

By induction a^_^ < min(a^~^, a^"*"^), so by the recurrence IT^ it 
follows that a^+i > max(a^~^, a^"*"^). (a) and (c) follow immediately. 

For the base case, we note that Oq = 1, while a\ = a^^ = c > 1, 
and so the result for Un follows. Notice that the above proof does not 
work for T„, since the base case fails. Indeed, if is the coefficient of 
in Tn{{c/2){x + 1/x)), then 6q = 1, while b\ = c/2, not necessarily 
bigger than one. However, we can use the result for Un, together with 
the observation (which follows easily from the addition formula for sin) 
that 

(13) Ux) = ^^(^^ 

Eq. (Uni) imphes that 6^ = - a^~^ > 0, by (c) above. □ 

The proof above goes through almost verbatim to show: 

Theorem 2.2. Let c > 1. Then all the coefficients of i?„ are non- 
negative. The same is true of Sn in place of Rn 

To complete the picture, we note that: 

Theorem 2.3. 

i?n(l;x) = i(x" + i^ 

Proof. Let x = expiO. Then l/2(x + 1/x) = cos9, and -R„(l;x) = 
T„(l/2(x + 1/x)) = cosne = l/2(x" + 1/x"). □ 

Remark 2.4. For c < —1 it is true that all the coefficients of Rn{c; .) 
and Sn{c; .) have the same sign, but the sign is (—1)"". For \c\ < 1, the 
result is completely false. For c imaginary, the result is true. I am not 
sure what happens for general complex c. 

By the formula (|TT| . we can write 
(14) 
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Noting that 

we obtain the expansion 
(16) 

n [2] 
Rn{c-x) = C^ E 

k=—n m=0 

where it is understood that (^) is if 6 < 0, or 6 > a, or 6 ^ Z. A 
similar formula for Sn can be written down by using Eq. (fT^. 



n 



n — m 



n — m 
m 



n — 2m 
{n-2m-k)/2j' 



3. Limiting distribution of coefficients 

While the formula ()16|) is completely explicit, and a similar (though 
more cumbersome) expression could be obtained for -R„(c; Xi, . . . , Xk), 
for many purposes it is more useful to have a limiting distribution for- 
mula as given by Theorem 13.11 below. To set up the framework, we 
note that since all the coefficients of Rn{c; Xi, . . . , Xk) are non-negative 
(according to Theorem 12. 2p . they can be thought of defining a proba- 
bility distribution on the integer lattice Z'^, defined by p{li, . . . ,lk) = 
[x^iX2 ■ ■ ■ x''i^]R{c] Xi, . . . , Xk)/R{c] 1, . . . , 1) (where the square brackets 
mean that we are extracting the coefficients of the bracketed mono- 
mial). Call the resulting probability distribution P„(c; z), where z now 
denotes a /c- dimensional vector. 

Theorem 3.1. With notation as above, when c> 1, the probability dis- 
tributions Vn{c; z/i/n) converge to a normal distribution on H'' , whose 
mean is 0, and whose covariance matrix C is diagonal, with entries 



2 



c 



^ k 



1 + 



, 1 \ 1/2' 
c + 1 ^ ' 



To prove Theorem 13. II we will use the method of characteristic func- 
tions (Fourier transforms), and more specifically at first the Continuity 
Theorem ( [Feller 1 1| Chapter XV.3, Theorem 2]), 

Theorem 3.2. In order that a sequence {Fn} of probability distribu- 
tions converges properly to a probability distribution F, it is necessary 
and sufficient that the sequence {4>n} of their characteristic functions 
converges pointwise to a limit (p, and that (j) is continuous in some 
neighborhood of the origin. 
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In this case is the characteristic function of F. (Hence is con- 
tinuous everywhere and the convergence (pn ^ (p is uniform on compact 
sets). 

The characteristic function 0„ of Vn{c] z) is simply 

Rn{c; exp{i9i), . . . ,exp{i9k))/ R{c] 1, . . . , 1). 
By definition of Rn, 

Rn{c;exp{i9i), . . .,exp{i9k)) = T„ Y!'j=i cos 9j^ , 

/?„(c;l,...,l)) =T„(fE'=iCos0) =r„(c). 
We now use the form of Eq. (jZj) : 



setting 

and 
we get 



u 



e 



Vcos^, 

j=i 



• • • , 9k), 



(17) 0. 



6> 



in J 2Tn{c) 

Notice, however, that for c > 1, the ratio of the second term in braces to 
the first is exponentially small as n — > oo, since the first term grows like 
(c + a/c^ — 1)", while the second as (c — ^ — 1)" (since cos ^ — > 1). 

Since, for the same reason, 2T„(c) = (c + \/(? — !)"[! + o(l)], we can 
write: 







n 



+ 0(1). 



Substituting the Taylor expansions for the cosine terms (hidden in u 
for typesetting reasons), we get: 

(18) u = k + ^{0,e) + o{l/n), 

In 
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SO 

(19) 



-u 



c + 



k 2kn 
A similar computation gives 



{e,e) + o{i/n). 



(20) 



Substituting the last expansion into the square root, we see that 



1 + ^ 



2n (c2-l)fc 

Adding Eq. (fTT?|l and collecting terms, get 
cu + V c^u^ — k"^ 
k{c+V^f^) 
1 



(c2-l)fc 



0.0)+o{l 



o{-) 



2nk 



1 



(C2 - 1)1/2 



{6,0) +oi-). 

n 



Performing some further simplifications, we see that 



e 



n 



exp 



U^ce\+o{i) 



where C is the covariance matrix described in the statement of Theorem 
13. H and Theorem 13.11 follows immediately. 
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